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Abstract 

For a given pseudo-Hermitian Hamiltonian of the standard form: H = p 2 /2m + v(x), 
we reduce the problem of finding the most general (pseudo-)metric operator r/ satisfying 
H< = r}¥Lr)~ x to the solution of a differential equation. If the configuration space is R, 
this is a Klein-Gordon equation with a nonconstant mass term. We obtain a general series 
solution of this equation that involves a pair of arbitrary functions. These characterize 
the arbitrariness in the choice of 77. We apply our general results to calculate 77 for the 
"PT-symmetric square well, an imaginary scattering potential, and a class of imaginary 
delta-function potentials. For the first two systems, our method reproduces the known 
results in a straightforward and extremely efficient manner. For all these systems we 
obtain the most general rj up to second order terms in the coupling constants. 
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1 Introduction 



The key aspect of the General Theory of Relativity (GR) that distinguishes it from other 
well-established physical theories is that in GR the very geometry of the spacetime, which is 
the arena of physical reality in classical physics, is itself a dynamical quantity. In contrast in 
Quantum Mechanics (QM) the geometry of the Hilbert space, which plays a similar role as the 
spacetime does in classical physics, is an absolute entity. This is to some extent dictated by the 
well-known mathematical fact that all separable Hilbert spaces are unitary-equivalent. For the 
past 75 years or so, this equivalence has been used to justify the absolutism associated with the 
convention of fixing the (inner product of the) Hilbert space from the outset. This is actually 
quite surprising, for the existence of an equivalence relation in a theory is clearly an evidence 
of the presence of a freedom in its formulation. In the case of QM, this is the freedom to choose 
the inner product of the Hilbert space, a freedom that has been left unused till recently, PIE]- 
In [2|, we have investigated the consequences of promoting the inner product of the Hilbert 
space into a degree of freedom. This revealed certain similarities between QM and GR and 
led to some interesting observations such as a direct link between geometric phases and the 
geometry of the Hilbert space and a new root to a certain nonlinear generalization of QM. In the 
present paper, we derive and examine a differential equation that includes among its solutions 
all possible choices of the inner product for a given physical system. This is the quantum 
mechanical analogue of Einstein's field equation. 1 For a system having R as its configuration 
space we obtain a series solution of this equation that involves two functional degrees of freedom. 
These signify the arbitrariness in the choice of the (pseudo-) metric operator. Our approach 
allows a more direct way of addressing some of the basic practical problems arising in the 
application of quasi- and pseudo-Hermitian quantum mechanics, |TJ |3J IH HE HI M M UH U21 UH] • 
In particular it provides an extremely powerful technical tool for the perturbative calculation 
of the (pseudo-) metric operators for various toy models. 

2 Differential Representation of Pseudo-Hermiticity 

Consider a physical system described by a separable Hilbert space 7i and a pseudo-Hermitian 
Hamiltonian operator H : 7i —>■ 7i. Let M. denote the set of all linear invertible Hermitian 
operators 7] : 7i — > 7i, then by definition |TJ] the pseudo-Hermiticity of H means that Ai H : = 
{i] G M.\H^ = r/Hi]^ 1 } is a nonempty subset of M.. The elements of M. are called pseudo- 
metric operators, for they may be used to define a pseudo-inner product (a nondegenerate 
sesquilinear form [T3]) (-|-) := (-\r}-) on H, where (-|-) denotes the defining inner product of 

lr The same way Einstein's equation does not generally restrict the metric tensor to have a particular signature, 
the above-mentioned equation does not restrict its solutions to correspond to positive-definite metric operators. 



2 



H 2 

The Hamiltonian H is Hermitian with respect to (-|-) for all 77 G M. H , |TT>^ I14j. If A4 H 
includes a positive-definite element, i.e., a metric operator, then H is Hermitian with 
respect to the positive-definite inner product (■\-) v+ - This implies that H is diagonalizable and 
has a real spectrum. The converse of this statement holds true at least for the case that the 
spectrum of H is discrete, i.e., if H is diagonalizable and has a real spectrum then Ai H includes 
a positive-definite element and equivalently H is Hermitian with respect to a positive-definite 
inner product, fTJEI]. Furthermore, in this case one can show that H is necessarily quasi- 
Hermitian, i.e., it may be mapped to a Hermitian Hamiltonian h : Ti — > Ti via a similarity 
transformation, H = p~ x hp, ^3^1- This and only this class of Hamiltonians are capable of 
supporting a unitary time-evolution in an associated physical Hilbert space. The latter is defined 
by endowing H with the inner product (-|-)^ and will be denoted by Tt v , [T§1 IH 151 17| 15]. 

An important fact about this construction is that 77+ is not unique. Different choices for 
i] + yield kinematically distinct quantum systems that nevertheless share the same dynamical 
structure. The quantum mechanical analogue of the principle of general covariance of GR is 
the physical (unitary-) equivalence of quantum systems , H), 0. The metric operators r/ + 
and more generally pseudo-metric operators rj are linked to and consequently determined by 
the Hamiltonian H via the pseudo-Hermiticity condition 

— r]Hr]~ 1 . (1) 

The same way Einstein's field equation links the metric tensor to the energy-momentum tensor, 
Q links the pseudo-metric operator to the Hamiltonian. The resemblance may be made more 
pronounced for a Hamiltonian of the standard form, 

H = ^ + v(x), (2) 

that acts in L 2 (R"). Applying both sides of ((TJ) on 77, substituting 0, and representing the 
resulting equation in the af-basis, we find 

9m \ 

- v ' + V S + -^K^r - v(y)}) v&y) = 0, (3) 

where V 2 := Y^j=i d 2 /du 2 for u = x,y and r](x,y) := (x\r]\y). For n = 1, this is a Klein-Gordon 
equation with a variable mass term, 

[-d 2 x + d 2 + /j 2 (x,y)]r ] (x,y)=0, fi 2 (x, y) := — [v(x)* - v(y)}. (4) 

According to (j2J) if r](x, y) is a solution, then so is i](y,x)*. The pseudo-metric operators 
rj G M. H correspond to solutions that satisfy 

v&y)* = v(y,x)- ( 5 ) 



^Strictly speaking, (-|-) is a nondegenerate sesquilinear form defined on the domain of rj. 
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Note that even for non-pseudo-Hermitian Hamiltonians of the form (J2J), Eq. (jSJ) admit solutions. 
However, these solutions fail to satisfy either the Hermiticity requirement (J^J) or the invertibility 
condition: 

/ d n y rj(x, y) ip(y) = implies tp = 0. (6) 
Jm™ 

If M. H happens to include positive-definite elements rj + , then these elements correspond to 
the solutions rj + (x, y) of ((3J) that in addition to (jSJ) satisfy 

d n x I d n y^{x)*r] + (x,y)ifj{y)>0 for ip ^ 0. (7) 

Jm™ 

The fact that for a diagonalizable Hamiltonian with a real and discrete spectrum such solu- 
tions exist is a consequence of the spectral theorems given in [El HHj. Indeed if {ip n , 0„} is a 
biorthonormal system associated with H, i.e., Htp n = E n ip ni H^(p n = E n (p n , (i/j n \4> m ) = S mn , 
then 

V+(%, V) = Y1 ^{x)<p n {yY (8) 

n 

is a solution of that satisfies both (j5J) and (J7J). As shown in (201 121]; in this case the 
most general (positive-definite) metric operator has the form A^i] + A where A is invertible and 
commutes with H. The latter corresponds to a solution of Q that is of the form 

V+&V)=Y,[ ^ [ d n vA(u,x)*cP n (u)MtYAv,y), (9) 
n Jm™ Jm™ 

where A (x, y) satisfies 

(-v* + vl + 2 -^[v(x) - v(jfi\) A&y) = o, (io) 



and 



/ d n y A(x,y)ip(y) = implies ^ = 0. (IT 

Jm.™ 



For a real- valued potential v, 5(x — y) is a solution of Q. It corresponds to the choice r\ = I, 
where / is the identity operator acting in 7i = L 2 (R). This is consistent with the Hermiticity of 
H . For a PT-symmetric potential v that satisfies v(— x) = v(x)*, 5(x + y) is a solution of (jSJ- 
This is a manifestation of P-pseudo-Hermiticity of the Hamiltonian [T3|, for (x\V\y) = 5(x + y). 

For n > 1, (jSJ) is an ultra-hyperbolic equation with quite peculiar properties [512] . We will 
therefore, focus our attention on the case n = 1. Our main purpose is to obtain the general 
solution of (j2J) without having to resort to its spectral decomposition flHJ). This is mainly because 
of the difficulties with summing the series in (jHJ) or evaluating the integrals that replace the 
latter whenever the spectrum becomes continuous [7j. 

In [0] , we have pursued a similar approach to construct the most general 7] + for the imaginary 
cubic potential v = iex 3 in some low orders of perturbation theory. The approach of [0] applies 
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to any (preferably imaginary) potential with a real spectrum. It yields an infinite system of 
iteratively decoupled partial differential equations whose solution provides the contributions to 
r] + in various orders of the perturbation theory. Although these equations have the same struc- 
ture, at each order one must compute their non-homogeneous term and solve them separately. 
In contrast, in the present paper, we obtain a single differential equation satisfied by 77, namely 
(j3J, that applies for an arbitrary potential v(x) rendering the Hamiltonian pseudo-Hermitian. 
An important advantage of the approach of the present paper over that of jU] is that in view of 
the simple structure of (j3J), we are able to offer (as discussed in Section 3) a general scheme for 
constructing a series solution of this equation. This solution involves two arbitrary functions 
that provide an explicit characterization of the arbitrariness in the choice of 77. 

Another recent application of the powerful machinery of differential equations to compute 
(pseudo-)metric operators is due to Scholtz and Geyer [TU]. These authors obtain a phase-space 
representation of the (pseudo-) metric operators r\. They use the Moyal product techniques to 
deal with the difficult factor-ordering problems that arise in this representation. The following 
are the main differences between the method of ^0] and the one presented in the present paper. 

• The method of (THj leads to an equation for r\ that is a differential equation [TT] provided 
that v (x) is a polynomial potential. Even for a polynomial potential the general character 
and in particular the order of this differential equation depends on the structure of v(x) 
and its degree. In contrast, in the present paper we offer a universal differential equation, 
namely (0J), that applies for polynomial as well as non-polynomial potentials, and has the 
same simple structure for all potentials. It is this appealing property that allows us to 
treat the well-known toy models of Section 4. The application of the method of ^H] to 
these models yields pseudo-differential equations (differential equations of infinite order) 
whose solution is extremely difficult if not impossible. 

• Suppose v is a polynomial potential, so that the method of JU| yields a differential 
equation, and suppose that one is able to solve this equation. Then one obtains an 
explicit expression for rj in terms of the operators x and p which involves a number of 
arbitrary functions. The condition that r] be Hermitian must be imposed to fix some of 
these functions. This is done by adopting a set of appropriate boundary conditions [XO j . 3 
In contrast, our method yields an expression for rj(x, y) that satisfies the Hermiticity 
condition i](x,y)* = r](y,x) manifestly and specifies a unique Hermitian rj according to 



It achieves this without making use of the Moyal product or having to select certain 

boundary conditions that ensures the Hermiticity of rj. Its successful application, however, 

3 The author is unaware of a systematic method of selecting the boundary conditions that achieve this 
purpose. 
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does not yield an explicit expression for 77 in terms of x and p. As explained in [9 , the 
latter may be obtained by Fourier transforming t](x,y) = (x\rj\y) over y to obtain (x\r]\p) 
and arranging the terms in the expression for y/2nh e~ lxp ^ h (x\r]\p) in such a way that x's 
are placed to the left of p's. This is how the issue of ordering of factors is addressed in 
this construction. 

A common feature of both methods is that solving the associated differential equations yields 
generally non-positive-definite pseudo-metric operators. The (positive-definite) metric opera- 
tors r] + , if they exist, correspond to certain special solutions that are to be identified using 
different means. 4 



3 Series Expansion for rj(x, y) 

We begin our analysis by expressing (jlj) in the form 

{-dl + d 2 y ) V (x, y) = f(x, y) - f(y, x)* (12) 

where 

2 772 

f(x, y) := v(y)r)(x, y). (13) 
We note that for a Hermitian rj, (|12j) is equivalent to 

(-d 2 x + dl) X (x,y) = f(x,y). (15) 

Next, we recall that the general solution of the wave equation (— + d 2 ) u(x, y) = is given 
by 

u(x,y) = u + (x-y) + u_(x + y), (16) 

where u± : ffi. — > C are a pair of arbitrary twice-differentiable functions (or distributions). 
Consequently, the general solution of (fT2jl has the form 

7](x, y) = u + (x - y) + u_(x + y) + x P (%, y) + Xp(v, x )\ ( 17 ) 

where u± satisfy u±(x)* = -u ± (=pr) and Xp{ x -> y) is a particular solution of (|15p. The latter is a 
non- homogeneous wave equation in 1+1 dimensions. It admits a particular solution which in 
view of (fT3*j) takes the form 



m ry px+y-r 

X P { x ,y) = T2 dr ds v(r)r)(s,r). (18) 

./ J x — y-\-r 



4 The construction of 77 given here may be supplemented with the procedure proposed in |1(J| for selecting 
the positive-definite metric operators rj + among 77's. This is expected to be a difficult task in practice, and we 
will not pursue it here. We suffice to point out that given r)(x, y) we can obtain an expression for 77 in terms of 
x and p as outlined in P| . This allows for making direct contact with the approaches of and |1 1 j . 
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Combining (JTfij) - (fTHj) . we find 

r](x, y) = u(x, y) + Kxi(x, y), 
where K is the integral operator defined by 



(19) 



m 



)Cr}(x,y) := — 



m 
ft? 



dr 



dr 



x+y—r 

x—y+r 
x+y—r 

x—y+r 



ds v(r) r](s, r) + I dr 



ds v(r) 7](s, r) + ds 



x+y—r 

■x+y+r 
x+y—s 

x+y+s 



ds v*(r) t](s, r)* 



dr v(s)* rj(s, r) 



• (20) 



In view of the analogy with the derivation of the Lippmann-Schwinger equation [23], it is not 
difficult to see that (fTUj) admits the following general series solution 



V(x, y) = [I - K] 1 u(x, y) = /C £ u(x, y). 



(21) 



£=0 



Clearly, rj is determined in terms of the arbitrary functions u±. 

For v = 0, i.e., a free particle, T](x,y) = u(x,y). As shown in 0, this is equivalent to 

77 = L(p) + K(p)V, 



(22) 



where L(p)' = L(p) and K(py = VK(p)V, equivalently L and K are respectively real-valued 
and VT- invariant 5 functions. 6 They are related to the Fourier transform 7 u± of u± according 
to 

L{p) = u+H), Kip) = V2^u4-l). (23) 

For a real-valued potential the ordinary choice for the metric operator that yields the L 2 - 
inner product, i.e., rj = I, corresponds to setting 

(24) 



ni 



u{x,y) = Six - y) - — J dr v(r). 

To see this, we first calculate KSix — y) for an arbitrary (possibly complex-valued) potential v. 
Using the well-known properties of the step function: 



9{x) 



for x < 
2 for x = 

1 for x > 0, 



(25) 



5 This means K(r)* = K(-r) for all rel. 

6 They may be further restricted to constants if one postulates the nonexistence of a hidden length scale for 
the problem. See [S] for details. 

7 In our convention, the Fourier transform of a function tp is given by (p{k) :— (27r) -1 / 2 dx e~ lkx <pix). 
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we then find 

( x-\-y x-\-y 

I 2 r 2 

/ dr dt[v(r)] + isign(y — x) / dr Qf[u(r)] ) , (26) 

where dt[v] and Q[v] respectively stand for the real and imaginary parts of v, and 

{—1 for x < 
for x = 
1 for x > 0. 

If f is a real potential, = i> and ^s[v] = 0. In this case (|26)) together with (|19|) and 
77(2; — y) — 5(x — y) yield IJ23J) . 

For a purely imaginary potential, 3?[w] = 0, v — i$s[v], and (J2fijl takes the following form. 

2 

dr f (r). (27) 

4 Applications 

4.1 PT- Symmetric Square Well 



PT-symmetric square well potential, 



— i( sign(x) for \x\ < 



L 



v(x) := S b W „ I (28) 

j oo tor \x\ > ij, 

with £ G R and L G M + , defines one of the best-known exactly solvable toy models that captures 
the generic properties of pseudo-Hermitian quantum systems [2H • A thorough investiga- 
tion of the physical content of this model is conducted in jlj where a particular perturbative 
calculation of a metric operator and the corresponding physical observables, localized states, 
probability density, and the underlying classical Hamiltonian is performed. This calculation 
makes use of the fact that the non-Hermiticity effects in this model diminish for energy states 
with larger spectral label N. More specifically it is (/N 2 that plays the role of the perturbation 
parameter. 

More recently, Bender and Tan ^2] performed a more conventional perturbative calculation 
of a metric operator taking ( as the perturbation parameter. This is the metric operator 
r] + that is associated with the CPT-inner product (-,-) CPT , That is (•, -) CVT = (■\r] + -), 

|21j . Expressing t] + in its exponential form, t] + = e~® , and noting that Bender and Tan set 
h = 2m — L/tz — 1, take e = — ( for the coupling constant, and use "e(x)" for " sign(x)", we 
can summarize their principal result (Eq. (11) of J2]) as 

(x\Q\y) =: Q(x, y) = -j[x-y+ sign(x - y)(\x + y\- tt)] + 0(( 3 ), (29) 
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where 0(( n ) stands for terms of order n and higher in powers of (. In particular, in view of 
the identity x — y = \x — y\ sign(x — y), we have the following expression for the CPT-metric 
operator r] + . 

it 

V+( X ,V) = S(x - y) + —(\x - y\ + \x + y\ - it) sign(a; - y) + C(C 2 )- (30) 

The perturbative calculation of the metric operator using the method developed in the 
preceding section is quite straightforward. Inserting ()28j) in (|27jl and performing the trivial 
integral in the resulting equation, we find 

K,5{x-y) = \x + y\ sign(ar-y). (31) 

The most general metric operator rj that reduces to the identity operator in the Hermitian limit 
£ — > is obtained by setting 

u(x, y) = 5(x -y)+ ([w+(x -y)+w_{x + y)\ + 0{C, 2 ) (32) 

in (j2U), where w± : [—§,§] — > C are arbitrary functions satisfying w±(x)* = w±(^fx) and 
w±(±L) = 0. 8 This together with (|HT|) yield 



r](x,y) =S(x-y) + C 



w + {x - y) + w-{x + y) + ^ \x + y\ sign(x - y) 



+ 0(( 2 ). (33) 



Setting h = 2m = L/ir = 1 in this equation, we find that the CPT-metric operator (|30j) 
obtained by Bender and Tan is a particular example of the metric operators ()33|) . It 
corresponds to the choice w + (x) = i(\x\ — n) sign(x) and w_(x) = 0. 

We can calculate higher order terms in the expression for the metric operator using our 
iterative method. Each additional order will involve an arbitrary pair of functions that enter 
the expression for u in (|2*T|) . This calculation is not only completely general (as it yields the 
most general metric operator), but it is also much simpler to perform. This is mainly because 
unlike its alternatives j3J it avoids approximating or summing complicated series. 

4.2 An Imaginary Scattering Potential 

Consider the following variant of the PT-symmetric square well potential |27j . 

«( , , — sign(x) for \x\ < ^ 

v \ x ) '■= 7T [sign(x + f ) + sign(a; - f ) - 2 sign(a;)J = <^ ; (34) 

2 I for |x| > ^, 

where £ 6 R is a coupling constant and L G M + is a length scale. 

8 These conditions arise from the Hermiticity requirement on the metric operator and its spectral resolution 
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In [7] we established the reality of the spectrum of this potential and used the spectral 
method of [T31 HZl UHJ to obtain a perturbative expression for an associated metric operator r] + . 
This involved constructing an appropriate biorthonormal system for the model and performing 
highly tedious calculation of the integrals appearing in the spectral resolution of r] + . Indeed, 
this calculation could only be done after expanding all the relevant quantities in powers of £ and 
restricting to the first order terms. Although the results reported in |7J required performing 
extremely lengthy calculations partly done using Mathematica, the expression obtained for 
i] + (x,y) := (x\r) + \y) took a surprising simple form, namely 

r] + {x,y) = 5(x-y) + ^ (2L + 2\x + y\ - \x + y + L\ - \x + y - L\) sign(x-y) + C(C 2 ). (35) 

Here, we wish to use the scheme developed in the preceding section to construct the most 
general metric operator rj that reduces to the identity operator in the Hermitian limit ( — > 0. 
In order to do this first we insert ()34|) in (|27J) and perform the trivial integral in the resulting 
equation to obtain 

JC5(x — y) = -—^ (\x + y + L\ + \x + y — L\ — 2\x + y\) sign(y — x). (36) 

Substituting in (j2TJ) and using (J31j|) . we then find 

r)(x,y) = 5(x-y) + ([w+(x-y)+w-(x + y) + 

(2\x + y\-\x + y + L\-\x + y-L\) sign(x -y)] + 0(( 2 ), (37) 

where w± : R — > C are arbitrary functions satisfying w±(x) = w±(±x). Clearly, the positive- 
definite inner product (|33j) constructed in [7j corresponds to setting w + (x) = sign(x) and 
W-(x) = 0. 



4.3 Imaginary (5-Function Potentials 

Consider the potential 

v(x) = i(5(x — a), (38) 

where (, a G M. 9 It is not difficult to solve the time-independent Schrodinger equation for this 
potential and show that H = 2— + v(x) has a real continuous spectrum. 10 This in turn suggests 
that one can construct an associated pseudo-metric operator using the spectral method of 
[T7| 11%] . This construction is similar to the one offered in [7] for the potential • An explicit 

9 Clearly, we can choose the origin of the x-axis so that a = 0. We retain a for future use where we consider 
the multi-delta-function potentials. 

10 H is not PT-symmetric. But one may attempt to use the results of [21] to construct a generalized VT- 
operator (an anti-linear involution) that commutes with H . 



10 



calculation of r\ using this method is however quite involved. A much simpler construction that 
we will describe in the following is based on the method of Section 3. 
First, we substitute (J3~%j) in (|2Uj) to establish 

KF(x,y) = F{x,y) + J 7 (y,x)*, F{x, y) := — 6(y - a) / ds F{s,a), (39) 

^ J x—y+a 

where z := 2m(/h 2 and F(x,y) is a test function. If we choose a ^-independent u, the series 
expansion (|2T|) becomes a power series in the coupling constant z. For defmiteness we shall 
first choose u(x,y) = 6(x — y). Setting F(x,y) = S(x — y) in (j3U|) and using the properties of 
the step function we then find 

%z 

Ku(x, y) = — 0(x + y - 2a) sign(y - x) =: m(x, y). (40) 

Alternatively, we could directly use (j2Zj) to obtain (I4~0l . 

Next, we compute K, 2 u(x,y) by substituting u\ for F in (|3l?|). This yields 

^ 2 m(x, y) = — - a) + 0(y - a)] [(x + y- 2a)9(x + y-2a)-\x- y|]. (41) 

The higher order terms in (|2*T)l can be similarly calculated. Moreover, because of the simple 
form of (j3T?j) and (j4T)j) . we can actually obtain an upper bound on \K l u{x, y)\ and use it to find 
a lower bound on the radius of the convergence of the series (pTj) . 

First, we recall |2H] that if a function g : K — ► C is bounded on an interval [a,/?] by some 
M G M + , i.e., |^(r)| < M for all r G [a,0\, then | rfr g(r)\ < M(j3 - a). Now, let F(x,y) 
be a function such that \F(s, a) \ has an upper bound M F as s takes values between \x — y\ + a 
and x + y — a. Then according to (}39|) . 

|/CF(x,y)| < |z|(|x-a| + |y-a|)M F . (42) 

In view of (|4f)|) . for all i,!/Gl, |/Cu(x, y)| < |z|/2. To obtain an upper bound on \JC 2 u(x, y)\ 
we set F(x,y) = JCu(x,y) in (J4*2*|) which allows us to identify M F with |-z|/2 and yields for all 

x,y G R: 

z 2 

\JC 2 u(x,y)\ < — (\x-a\ + \y - a\). (43) 



We can directly verify this relation using 1)41)1 . Repeating the procedure that leads to 
we find for all t > 1 and all x, y G R: |/C £ w(a;, y)| < |^| (|x — a| + |y — a|) f_1 /2. This in turn 
implies, in view of the elementary comparison tests, that the series 1)21)1 converges (absolutely) 
for I^KIx — a| + |y — a\) < 1. Hence, it converges in an open disc in the x-y plane that is centered 
at (x — a,y — a) and has a radius g > [\[2 |-2|) _1 . 

In summary, for every given value of z, truncation of the series ()21)) yields a reliable approx- 
imation for 7] provided that we keep sufficiently large number of terms in the series and deal 
with wave functions ip(x) that decay sufficiently rapidly as \x\ — * oo. 11 



ii 



Our analysis only yields a lower bound on g. It does not imply that g is finite. 
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We can extend our treatment to a potential consisting of more than one delta-function: 

N 



J2(nS(x-a n ) } (44) 



V{X) = I 

n=l 

where iV G Z + and ( n , a n G R. An example is the PT-symmetric potentials [2H1 EH EH] 
corresponding to the cases that N is even and (i± +k = a^+k — ~ a k for all k — 1, 2, • • • y- 12 
For these multi-delta-function potentials the calculation of the first order term in z n : = 
2m( n /h 2 in the series expansion (J2~T]) reduces to the case N — 1 that we considered above. In 
view of (f47Ijl . 

. AT 

0(x, y) = - ^2 z n 9{x + y - 2a n ) sign(y - x). (45) 

n=l 

The most general rj that reduces to 77 = J in the Hermitian limit z n —> is obtained up to 
second order terms in z n by setting u(x, y) = S(x—y)+^2 n _ t z n [w n+ (x— y)+w n -(x+y)}+0(z 2 ) 
where w n ± : R — > C are arbitrary functions satisfying w n± (x)* = w^^x). This together with 
flUD and (j2U) yield 



77(x,y) = 5(x-y) + 



N 
n=l 



w n+ {x -y) + w n ^(x + y) + -0(x + y - 2a n ) sign(y - x) 



+ 0{zl) 



(46) 

We close this section by the following general remarks. As we observe in the study of the 
above toy models, the series solution pTf may be used to obtain a perturbative expansion for 
the pseudo-metric operator 77. In general depending on the details of the model under study 
one may or may not have access to a dimensionless perturbation parameter. Typical examples 
for which this occurs are the imaginary cubic potential and the single imaginary delta-function 
potential (JHHj) . In this explained in P , the truncation of the perturbative expansion of 

rj(x,y) generally yields a reliable result only within a sufficiently small region in the x-y plane. 
Furthermore, one expects that for sufficiently small values of the perturbation parameter (the 
coupling constant ( or £ ra in the above examples) the perturbative corrections to a positive- 
definite metric operator such as 77 = / leave this property intact. 



5 Concluding Remarks 

In this article, we have outlined a differential realization of the pseudo-Hermiticity condition 
that plays a central role in devising a unitary quantum theory based on quasi- Hermitian Hamil- 
tonians of the standard form. The integral kernel r](x,y) for the corresponding pseudo-metric 

12 Here we assume that ((„, a n ) are such that the Hamiltonian is pseudo-Hermitian. This is the generic case, 
for the values of (C, n ,a n ) that render the spectrum of the Hamiltonian non-real form a measure-zero subset of 
the set M. 2N of all possible values of a n ). 
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operators r\ satisfies a linear partial differential equation. For systems having K. as their config- 
uration space this is nothing but a particular variable-mass Klein-Gordon equation. We have 
obtained a general series solution for this equation and demonstrated its application in treating 
the PT-symmetric square well potential, an imaginary "PT-symmetric scattering potential, and 
a class of imaginary delta-function potentials. In particular for the former two potentials, the 
approach presented here is by far more practical than the alternative approaches that use the 
spectral resolution of the metric operator. Another advantage of the former approach is that it 
is capable of producing the most general pseudo-metric operator. In particular, imposing the 
positive-definiteness condition (J7J), it yields the general form of the metric operators. 

Our method is not only practically advantageous but also conceptually appealing. It furthers 
the analogy between Quantum Mechanics and General Relativity, for the differential pseudo- 
Hermiticity relation plays a similar role in Quantum Mechanics as the Einstein's field equation 
does in General Relativity. 13 Another valuable outcome of our method is a concrete charac- 
terization of the arbitrariness of the metric operator. Each choice of a metric operator defines 
a separate quantum system. One can pursue the prescription used in the so-called quasi- 
Hermitian quantum mechanics pQ to select an irreducible set of compatible quasi- Hermit ian 
operators O a and fix the metric operator 77+ (up to scale) through the requirement that O a 
be ?7 + -pseudo-Hermitian. Alternatively, one can follow the approach of the so-called pseudo- 
Hermitian quantum mechanics choose 77+ directly, and construct the Hilbert space and 
observables of the theory accordingly |U E] . 
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